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Abstract. A new class of high mass X-ray binaries (HMXBs) — supergiant fast X-ray transients (SFXTs) — are discovered 
by INTEGRAL, which are associated with OB supergiants and present long spin periods. There are observational evidences 
indicating that some accreting neutron stars in HMXBs display accretion reversals. It has been suggested that the inverted 
torque can lead to a very slow rotator. According to three characteristic radii in wind-fed accretion, we develop a retrograde 
accretion scenario and divide the accretion phase into three regimes, in order to interpret the formation of the long spin period 
of SFXTs. The accretion regime in some SFXT systems can be determined by their spin and orbital periods. 
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1. Introduction 

Supergiant fast X-ray transient (SFXT) is a new class of mas- 
sive X-ray binaries recognized by INTEGRAL observations 
(Sguera et al. 2005; Negueruela et al. 2006), which contains a 
magnetized neutron star accreting material from an OB super- 
giant companion star. They, showing short and sporadic out- 
bursts (Romano et al. 2007) and spanning 2-5 orders of mag- 
nitude in dynamic range (e.g. in't Zand 2005; Leyder et al. 
2007), represent the extreme case of X-ray variability in high 
mass X-ray binaries (HMXBs) with early-type supergiant com- 
panions (Masetti et al. 2006b; Pellizza et al. 2006). The quies- 
cent state is characterized by a low luminosity in the range of 
~ 10 31 - 10 33 ergs/s (e.g. in't Zand 2005; Leyder et al. 2007). 

Some SFXTs are believed to be X-ray pulsars with long 
spin periods, and their measured pulse periods are located in 
a wide range, e.g. from 4.7 s (AX J1841.0-0536/IGR J18410- 
0535, Bamba et al. 2001) to 1240 s (IGR J16418-4532, Walter 
et al. 2006; Sidoli 2011a). Generally speaking, a newly born 
neutron star is presumed to be a fast rotator with a spin period 
of a fraction of one second. If a neutron star is in a binary sys- 
tem, it may experience an ejector state (the conventional spin- 
powered energy-loss mechanism) and propeller state (the in- 
teraction between magnetosphere of the neutron star and the 
stellar wind of the companion) before accreting material from 
its optical companion (for a review see Bhattacharya & van 
den Heuvel 1991, and references therein). During the ejector 
and propeller regimes, the neutron star spins down and forms a 
slow rotator (van den Heuvel 1977; Davies, Fabian & Pringle 
1979; Davies & Pringle 1981) at the end of these two phases. 



Under certain assumptions (see Urpin, Konenkov & Geppert 
1998; Li & van den Heuvel 1999), the spin period can be lim- 
ited to a few hundred seconds (~ 500 s). If the system under- 
gos a spherical accretion and three epochs, i.e. ejector, super- 
sonic propeller and subsonic propeller epochs, the conditions 
under which a longer period can be achieved may be relaxed 
(Ikhsanov 2001 ; Ikhsanov 2007). 

Long-period pulsars are often identified as an accretor and 
their periods gradually decreases during the accretion in the 
prograde scenario. However, observations indicate that dy- 
namic torque reversals are present in some wind-fed persistent 
HMXBs (Bildsten et al. 1997; Nelson et al. 1997, and refer- 
ences therein). The spin-down torque due to the accretion rever- 
sals has the same order of magnitude as the spin-up torque, but 
opposite in sign (Nelson et al. 1997). Accretion torque reversals 
can occur in an inhomogeneous wind (Ruffer 1997, 1999). 

In this paper, we suggest that the long spin periods are as- 
sociated with a retrograde accretion scenario. In section 2, we 
review the theory of wind accretion in SGXBs/HMXBs and de- 
rive the spin evolution and the spin-down timescale in different 
accretion regimes. We apply our model to SFXTs in section 3. 
Section 4 contains a short summary. 

2. Retrograde wind accretion 

OB supergiants can produce powerful stellar wind driven by 
radiation, which is approximately isotropic, that is, the wind 
material leaves the surface uniformly in all directions, with a 
wind velocity of v w ~ 1000-2000 km s _I . The OB companions 
in SFXTs typically have the mass of M c ~ 30M Q , mass loss 
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rate of M w ~ 10~ 6 M© yr~ l , and luminosity of log(L c /L©) ~ 
5-6. Here, M© is the solar mass, and L© is the solar luminosity. 

The magnetized neutron star in SGXBs accretes matter 
from an OB companion via the intense stellar wind (Verbunt 
& van den Heuvel 1995). We consider the case of spherically 
symmetric stellar wind radiated from the OB companions and 
disordered and inhomogeneously captured by relativistic neu- 
tron stars, and the accreting neutron star can be either a pro- 
grade accretor or a retrograde accretor (Bildsten et al. 1997; 
Shakura et al. 2012). If the neutron star behaves like a retro- 
grade accretor, it may spin down and rotate with a very long 
period because of the inverted accretion torque. 

2.1. Characteristic radii 

In the theory of wind accretion by a magnetized rotating neu- 
tron star, three characteristic radii (Davidson & Ostriker 1973; 
Lamb, Pethick & Pines 1973; Stella, White & Rosner 1986) can 
be defined as the following. Only part of the wind material from 
the OB companion will be captured by the gravitational field of 
the neutron star. The radiated wind flowing towards the neutron 
star within a certain radius will be accreted, whereas the flow- 
ing material outside that radius will escape. This radius, called 
the accretion radius r acc (Bondi & Hoyle 1944), can be defined 
by noting that material will only be accreted if it has kinetic 
energy less than the potential energy in the gravitational field 
of the neutron star, that is, 



M 



(1) 



racc = 2GM/vi~W"— v-%, 

where M is the neutron star mass, and v w> g is a velocity in units 
of 1000 km s- 1 (v w , 8 = v w /(10 8 cm s" 1 )).' 

When the magnetic pressure of the neutron star can balance 
the ram pressure of the accreting material, a magnetosphere 
forms and the inflowing matter may stop at such a preferred 
radius — magnetospheric radius r mag (Eisner & Lamb 1977; 
Davies & Pringle 1981). It can be defined by 

B(r mag ) 2 



• p(.Fmag)v(.fmag) ■> 



(2) 



r \' magJ v V mag. 

where B(r mag ) ~ nlr^ ag is the magnetic field strength at r mag , 
and the p(r mag ) and v(r mag ) are wind density and wind veloc- 
ity at r mag , respectively, fi is the magnetic moment of the neu- 
tron star. In the case of r mag > r acc , the inflows cannot expe- 
rience significant gravitational field. So the wind speed at r mag 
is v(r mag ) = v w . Assuming a nonmagnetized spherically sym- 
metric wind (Eisner & Lamb 1977), the wind density p(r mag ) 
reads p(r mag ) « M w /4na 2 v w (Davies & Pringle 1981; Bozzo, 
Falanga & Stella 2008), where a is the orbital separation of the 
two components in the system and that a » r mag is assumed. 
Therefore, the magnetospheric radius is 



r mag ~ 10 10 AT 1/6 "" 1/6 " 1/3 " 1/3 



(3) 



where M w ,_ 6 = MJW+MJyi), // 33 = y"/(10 JJ G cm j ), 
and «io is related to the orbital separation as a ~ 10 12 flio cm 
~ IO^P^jqM^ 3 cm. P h jo is the orbital period P b in units of 
10 days, and Mjq is the companion mass in units of 30 so- 
lar masses. When r mag < f acc , the gravitational field of the 



neutron star dominates the falling of wind matter. The wind 
density and wind velocity at r mag can therefore be taken as 



P^Xmag) 



and v(r mag ) 



v(r mag )A7ir; nag 

accretion rate M depends on r a 
King & Raine 2002) 

M/M w ~ r 2 ee /(4a 2 ) ~ l(T 5 (^-) z v 

Then the magnetospheric radius is given by 

r „ 10 1( Y— r 5/7 Ar 2/ V /7 a 4/ V /7 cm 

rmag M© M w,8 10 ^33 U " 



-, respectively. The 



2GM 

y fmag 

and a according to (Frank, 



w,& a W ■ 



w,—6 w,8 1 



(4) 



(5) 



A corotation radius r cor can be defined by the spin angu- 
lar velocity (Q. s = 2n/P s ) of the neutron star being equal to 

the Keplerian angular velocity (Q.k = yjjr-) of matter being 

accreted, 



M 



M© 



,2/3 
.5,3 



cm 



(6) 



where P s ^ is the spin period of the neutron star P s in units 
of 10 3 s. The competition among the three characteristic radii 
leads to different accretion regimes, so do the accretion modes 
and the behavior of the neutron star. 

2.2. Pure propeller 

In a system in which the magnetospheric radius is larger than 
the accretion radius (r mag > r acc ), the stellar wind, flow- 
ing around the magnetosphere, may not experience signifi- 
cant gravitational field and interacts directly with the magne- 
tosphere. The neutron star behaves like a pure propeller and 
spins down via rotational energy dissipation resulting from the 
interaction between the magnetosphere and stellar winds, cor- 
responding to radio pulsars with magnetic inhibition (Stella, 
White & Rosner 1986; Bozzo, Falanga & Stella 2008). The 
spin-down torque can be expressed as (Lynden-Bell & Pringle 
1974; Wang 1981; Lipunov 1992) 

Tp, P ro — —KtH /r mag > (7) 

where K t < 1 is a dimensionless parameter of the order of unity 
(Davies & Pringle 1981; Lipunov 1992). The equation govern- 
ing the spin evolution can be written as (Davidson & Ostriker 
1973; Ghosh & Lamb 1978) 



„ d 1 
2n1—— = T , 
dtP, 



(8) 



where / is the moment of inertia of the neutron star, and T 
is the total torque imposed on the neutron star. Therefore, the 
spin-down rate reads 

P s , pur ~ lO-^p^P^a-^yj 2 , (9) 

where /4s is the moment of inertia in units of 10 45 g cm 2 . The 
spin-down timescale is given by 

T s , pur = P s /P s , pU r ~ 10 3 ^ 1 ^ 33 1 /45P;>10M; 1 _ / 6 2 vrv! 8 /2 - ( 10 ) 

The luminosity in this regime results from the shock forming 
near r ma „, 



r 2 H 

magr 



(11) 
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2.3. Retrograde propeller — quasi-spherical capture 2.4. Retrograde accretor — quasi-disk accretion 



If the magnetospheric radius is smaller than the accretion ra- 
dius, the stellar wind radiated from the companion can pene- 
trate through the accretion radius and halt at the magnetopause. 
We consider the regime of r mag > r cor in this part. The inflow in 
the region between r acc and r mag falls approximately in a spher- 
ical configuration. The neutron star rotates at a supersonic ve- 
locity and ejects some of that wind material out of the accretion 
radius, leading to the dissipation of some rotational energy and 
thus spinning down. Such a scenario corresponds to the pro- 
peller mechanism and imposes a spin-down torque on the neu- 
tron star. The torque resulting from the propeller mechanism 
reads (Lipunov 1992) 



T r , pro = -Mvl/Q s . 



(12) 



In addition, if the accumulation of accreted matter is faster 
than the ejection rate by the propeller mechanism, more and 
more wind material may deposit outside the magnetosphere 
(Stella, White & Rosner 1986). Because of the disorder of the 
stellar wind and the turbulent motion resulting from the su- 
personic rotation of the neutron star (Davies & Pringle 1981), 
the flowing accumulated material can be either prograde or 
retrograde around the magnetosphere. If the wind flows ret- 
rogradely, it will impose an inverted accretion torque on the 
neutron star and spin it down. Such an accretion torque can be 
written as (Lipunov 1992) 



TV 



acc — M^l s r ma g 



(13) 



Accordingly, the total torque imposed on the neutron star 
in this regime is 



T P = T r , pro + T nacc = -Mv 2 JCl s - MQ, s r, 



mag 



Substituting Eq. (TBI into Eq. (|S), we can obtain the spin evo- 
lutionary law in the regime of r acc > r mag > r cor , 



M 



Kpro ~ 10-^(^ w ,-6V;Xo 2 ^,3 + (^) 

Therefore the spin-down timescale reads 



9/7 



(15) 



10V45(M w ,_ 6 v-Wp2 + (__)' 



M 



,9/7 



Mr, 



M 5P <T 10/ V 20/ V /7 ) vr 

m w,-6 u W v w,S ^33 > J L 



(16) 



Both the interaction with magnetic fields by the propeller 
mechanism and inverted accretion contribute to the luminosity, 
which reads 



Lpro = Mvl + GMM/r, 



mag 



M 



• >9/7 -36/7 -18/7 -4/7-. -1 

M J,-6 v w,% a w ^33 )ergss . 



(17) 



In addition, the wind material halts at the magnetopause and 
forms a shock near r mng . Energy is also released through the 
shock, and the corresponding luminosity is given by Eq. (17) 
in Bozzo, Falanga & Stella (2008). 



When r acc > r mag and r cor > r mag , the disordered stellar wind 
penetrates into accretion radius and falls directly toward the 
magnetosphere. Owing to the relatively slow rotation of the 
neutron star, the rotating magnetosphere cannot eject the wind 
material. The inflows experience significant gravitational field 
and accumulate near the magnetosphere. The accumulated ma- 
terial, with high initial velocity, can form an extended adiabatic, 
tenuous and retrogradely rotating atmosphere around the mag- 
netosphere, imposing an inverted accretion torque and spinning 
down the neutron star. In addition, the reversal rotation between 
the adiabatic atmosphere and the magnetosphere may cause a 
flip-flop behavior of some clumps in the atmosphere, ejecting a 
part of material and leading to the loss of thermal energy. 

The spin-down torque resulting from the flip-flop ejection 
has the same form as Eq. (0 (Davies & Pringle 1981; Lipunov 
1992). The accretion spin-down torque imposed by the inverted 
atmosphere has the same expression as Eq. (1131 1, So the total 
spin-down torque imposed on the neutron star in the regime of 



and r,, 



T a = -k,-j MCl s rl ag 

f"mag 



(18) 



Substituting the above expression into Eq. ([8j, we get the spin 
evolution equation, 



M 



P ~ 10- 6 rV^/; 2/7 P 2 ( — ) l5 l 1 M 6 ' 7 y- 24/ V 12/7 
r S ,acc '45W33 '.^M. 10 



M 



+ ^/ 7 p„3(^)- i4/7 <i 6 v: 6 ^^ 7 ) - s- 1 ■ d9) 



(14) The spin-down timescale is 



M 



Ts , acc ~ io 3 / 45 (^p. s ,3(-^) I ^<I 6 v- 4 / 7 fl - 1 ^ 



+ ^ 7 (^)- 14/7 <I 6 « 7 )" 1 yr 

The corresponding luminosity is 



M 



Lace ~ 10 31 (^) 26/7 <I 6 v-^- 8/ V33 4/7 ergs s-i 



"M, 



(20) 



(21) 



3. Application to SFXTs 

Some SFXTs containing sporadically accreting neutron stars 
are believed to be X-ray pulsars and display very different pulse 
periods, from a few seconds (e.g. IGR J16479-4514, P s = 2.14 
s, Sguera et al. 2008; IGR J18483-0311, P s = 21 s, Sguera, 
Hill & Bird et al. 2007; IGR J17544-2619, P s = 71.49s, Drave 
et al. 2012; AX J1749. 1-2733, P s = 132 s, Karasev et al. 2007, 
Karasev et al. 2008; IGR Jl 1 1215-5952, P s = 187 s, Swank 
et al. 2007; IGR J16465-4507, P s = 228 s, Lutovinov et al. 

2005) to a few thousand seconds (e.g. IGR J16418-4532, P s = 
1240 s, Walter et al. 2006, Sidoli et al. 2012; XTE J1739-302, 
reported for periodicity in the 1500-2000 s range, Sguera et al. 

2006) . Their orbital periods also spans over a very wide range, 
from 1.4 days to a few hundred days (e.g. Sidoli 2011b, and 
references therein). 
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According to the definition of the three aforementioned 
characteristic radii (see Eq. (fl}, Eq. (|5j and Eq. (O), they de- 
pend on some physical properties of the system, including the 
neutron star mass and magnetic field strength, stellar wind ve- 
locity, orbital period of the binary system, the mass loss rate 
of stellar wind, the mass of the companion and the neutron 
star spin period. In this work we take the neutron star mass 
as IAMq and magnetic field as 10 12 G, the mass of OB su- 
pergiant as 3OM , the mass loss rate and velocity of the stellar 
wind as 1O~ 6 M yr~' and 1000 km/s, respectively. Then the 
accretion radius is about 10 10 cm, and the corotation radius de- 
pends only on the spin period as r cor ~ P 2 J^- The magneto- 
spheric radius depends also only on the orbital period accord- 
ing to r mag ~ P 2 ^ in the region of r mag > r acc and r mag ~ P*^ 
when r mag < r acc . Therefore, we can determine the accretion 
regime in which a certain SFXT is located, according to the 
spin period and orbital period (see Table 1). 

IGR J18483-0311 and IGR J17544-2619 display not-too- 
slow rotation and relatively short orbital periods (for IGR 
J18483-0311, P s = 21 s, Sguera et al. 2008, P orh = 18 days, 
Jain et al. 2009; for IGR J17544-2619, P s = 71.49 s, Drave et 
al. 2012, P orb = 4.9 days, Clark et al. 2010). From the expres- 
sions of r mag and r cor (see Eq. (0 and Eq. (0), one can see that 
the relation among these three radii is r acc > r mag > r cor . So 
IGR J 1848 3 -03 1 1 -like sources can behave like a propeller in 
the retrograde scenario. 

Very long spin period SFXTs with narrow orbits, e.g. IGR 
J16418-4532 (P s = 1240 s, Walter et al. 2006, Sidoli et al. 
2011a; P or b = 3.74 days, Levine et al. 2011), has a typically 
large corotation radius. Moreover, the magnetospheric radius 
is one order of magnitude smaller than the accretion radius. So 
this class of systems can be an accretor with retrograde accre- 
tion. In addition, IGR J16479-4514, with P s = 2.14 s (Sguera 
et al. 2008) and P or i, = 3.3 days (Jain et al. 2009), also presents 
a mildly larger corotation radius. Therefore, it can also behave 
like an retrograde accretor. 

There are several systems rotating slowly with spin peri- 
ods of a few hundred seconds and relatively wide orbital sep- 
aration, such as IGR Jl 11215-5952 (P s = 187 s, Swank et 
al. 2007; P orb = 165 days, Sidoli et al. 2007), AX J1749.1- 
2733 (P s = 132 s, Karasev et al. 2007, Karasev et al. 2008; 
P orh = 185.5 days, Zurita & Chaty 2008), and IGR J16465- 
4507 (P s = 228 s, Lutovinov et al. 2005; P orh = 30.32 days, 
Clark et al. 2010). As listed in Table 1, the relations among the 
three aforementioned characteristic radii of these sources are 
rmag > face > r cor- Therefore, these three sources can be in the 
pure propeller regime. 

The outbursts of XTE J1739-302 observed by several ob- 
servatories show a characteristic timescale for variability of the 
order of 1500-2000 s (Sguera et al. 2006). This source has an 
orbital period of 51.47 days (Drave et al. 2010). Both the coro- 
tation radius and magnetospheric radius are mildly larger than 
the accretion radius (see Table 1 for details). However, due to 
the uncertainty in the spin period, it is difficult to determine its 
exact status. More accurate knowledge about the properties of 
the system is needed to determine whether the system is a pure 
propeller or a retrograde propeller. 



It can be seen that the three characteristic radii in several 
sources have the same order of magnitude (see Table 1 for de- 
tails). IGR J16418-4532 shows a dynamical range of two orders 
of magnitude and quasi-periodic flares (Sidoli 2011b). These 
features were taken to suggest that its X-ray emission is driven 
by a transitional accretion regime in a prograde accretion sce- 
nario (Sidoli 2011b). In this sense, we claim that this class of 
SFXTs frequently experience a transitional state, between the 
propeller regime and the retrograde accretor phase. 

The above discrimination depends strongly on the assump- 
tion that these systems have the same properties except for spin 
and orbital periods. The particularities for each system should 
be investigated for further detailed study. 

4. Summary and Discussion 

Neutron stars in a wind-fed accreting system with disordered 
and inhomogeneous stellar winds can be either a prograde or a 
retrograde accretor. An accretor may exhibit different behaviors 
in different accretion regimes, depending on the competition 
among the three characteristic radii. 

When the magnetospheric radius is located outside the ac- 
cretion radius, the wind cannot experience significant gravi- 
tational focusing and forms a bow shock near the magneto- 
sphere. The interaction between the stellar wind and the mag- 
netosphere contributes to the loss of rotational energy and the 
spin-down of neutron star. This is the pure propeller regime. 

If the accretion radius is larger than the magnetospheric ra- 
dius, the wind flowing from OB companions with a high veloc- 
ity interacts with the magnetosphere directly and forms a shock 
at rmag (Harding & Leventhal 1992). In such a case, when r mag 
is larger than the r cor , the neutron star can behave like a ret- 
rograde propeller. The shock, formed near the magnetosphere, 
flip-flops from side to side and makes the uniform wind disor- 
dered and inhomogeneous. Therefore, the neutron star accretes 
matter with a relatively high angular velocity, but with opposite 
signs (Livio 1992), consistent with the retrograde scenario. 

If both the accretion radius and corotation radius are 
larger than the magnetospheric radius, the inflowing matter can 
fall directly and halt near the magnetopause. Because of the 
high falling velocity, there is no effective cooling mechanism. 
Moreover, the dissipation of the rotational energy also produces 
additional heat. The spherically symmetric accreting material 
forms a quasi-disk, adiabatic and tenuous atmosphere around 
the neutron star, which can retrogradely rotate with the magne- 
tosphere and spin down the neutron star. 
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